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Abstract 

We study in this article the mathematical properties of a class of orbital-free kinetic energy 
functionals. We prove that these models are linearly stable but nonlinearly unstable, in the sense 
that the corresponding kinetic energy functionals are not bounded from below. As a matter of 
illustration, we provide an example of an electronic density of simple shape the kinetic energy of 
which is negative. 

1 Introduction 

Kohn-Sham models have brought a considerable breakthrough in atomic-scale simulation of materials in 
condensed phase. However, the use of the Kohn-Sham kinetic energy is problematic when the number of 
electrons by unit cell exceeds a few hundreds (for computational means available to date) . Some authors 
therefore proposed to approximate the Kohn-Sham kinetic energy functional in order to get rid of both 
the orbital and /c-point dependencies. Their approach consists in improving the Thomas-Fermi model, 
for which the kinetic energy functional reads 



TtfW = Ctf / p'^' (1.1) 

3 2/3 

where Q is the simulation unit cell, p a given density, and Ctf = l"^^^) Thomas-Fermi constant, 

by adding some correction terms. The functionals under consideration in this article are referred to as 
density-independent in the literature jH]. They formally read 

T„,/5[p] - Ctf / p'^' + l f \^Vp\' + Ctf [ ( [ w^,f,{ko[p],x - y) p{y f dy) dx (1.2) 

JQ ^ Jq Jq \Jwt> J 
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1 /3 

where a and /3 are positive real numbers such that a + f3 = 5/3, where fco[p] ~ (S^'^p) is the Fermi 
wavenumber associated with the average density p = --— f p (here and below \Q\ denotes the volume 

\Q\ Jq 

of the unit cell Q), and where Wa,f3 is some Green kernel. We will denote respectively by Ttf, Tw and 
Tc the Thomas-Fermi, von Weizsacker and convolution term in lll.2|l : 



TtfM = Ctf / p'/', (1.3) 

JQ 

Tw[p] IVVpI' (1.4) 

Tc[p] = Ctf l^pixT (^l^^w^Ako[p],x-y)p{yy'dy^ dx. (1.5) 

Note that Wa,0 is a function of two variables. The first one, denoted by fc^, is a real number which has the 
dimension of a wavenumber. The second one is the convolution variable; it is a vector of which has the 
dimension of a position vector. Energy functionals of this type were introduced by Wang and Teter \^ 
(with a — P — 5/6), and further generalized by several authors j^lCIIHI- For a given pair (a, /S), the Green 
kernel Wa^fs is completely determined by the requirement that the kinetic energy functional T^jj must 
be compatible with the Lindhard perturbation theory (see e.g. ^). This compatibility condition has 
been written as early as in 1964, in the article by Hohenberg and Kohn founding the Density-Functional 
Theory p]. Imposing that Ta^p must be compatible with the Lindhard theory leads to the relation 

<''(^-^)= 9^^(11) 

where Wa,/? denotes the Fourier transform of Wa,f3{kF,x) with respect to the convolution variable x and 
where for all ?/ G K+, 

■ 1 + r/ 



2,1 f ~ 1. 



It is important to note that the normalization convention entering in the definition of the Fourier transform 
used above is the following: for all / G L^( 



The purpose of this article is to analyze the mathematical properties of the kinetic energy functionals 
of the form ljl.2ll . The main results are presented in Section [2 We prove that these models are linearly 
stable but nonlinearly unstable, in the sense that the corresponding kinetic energy functionals are not 
bounded from below. As a matter of illustration, we provide an example of an electronic density of simple 
shape the kinetic energy of which is negative (all the numbers are in atomic units). Let us consider a 
cubical simulation cell Q —] — i/2,L/2p and the Q-periodic function PN.,ro,L with N > Q and tq > 
defined on Q by 

p^,„,,(x,y,z) = iv (lS'\-{^'+y'+^')/rl. 
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For rQ << L, one has (up to machine precision), 

PN.r„,L — N, 



CtF /3\'/' 5/3 



Ttf [PN,ro,L] - j ^° ' 

and 

5 

In the above expressions, Ttf, Tw and Tc are defined by lll.3|l . ljl.4p and ljl.5ll . and 

M7) = ^h^ E G(7|g|)e-I^l^ (1.6) 

ge(/iZ)3 

with 

5 V^' ro /3a/3\^/^ 1 

TT— and 7 



Notice that is a Riemann sum which approximates the integral 

'^(t)-4j / G(7k|)e^'^l^dg. 

For — 13, To = 0.5 and L — 4.906, and with a, /3 = — 6~~~' instance Tk^ipAr.ro.L] = —8.183 

(for the sake of comparison, the kinetic energy of the uniform electron gas of density p — N/L^ is 
Ta,p[p] ~ 8.573). The parameters = 13 and L = 4.906 correspond to an all electron calculation on 
Aluminium. 

2 Main results 

Let us first recall the definitions of the functional spaces under consideration below: for 1 < p < +oo, 
LP ^(M^) = |m : m3 ^ M measurable, J \u\p < +oo for all compact sets K cM.^ | , 

-^^(Q) = |u : Q ^ M measurable, J \u\p < +cx)| , 
H\M.^) ^ {ii e L'^iR^), Vii e (l2(r3))^} , 
H^M = {^^= v\q, t'eiL(R'), e {Ll^{R')f , t; g-periodic} , 
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-f^pcilQ) = i^lcriQ))' is the topological dual of Hp^j.{Q). 
The space U'iQ) is equiped with the norm = (^j and the space Hp^^.{Q) with the norm 



1/2 

Mm = ( / \u\'+ I \Wu\' 



The main mathematical properties of the orbital-free kinetic energy functional Ta.j3 are put together in 
the following two theorems. 

Theorem 2.1 Let us consider two positive real numbers a and (3 such that a + (3 — 5/3. Let us consider 
a Q -periodic potential V e L'^l^iM?) and the minimization problem 

In = inf I^T^Ap] + J^Vp, p > 0, ^eHl^^{Q), ^ p - ivj , (2.1) 

where N is the number of electrons per unit cell. Then, 

1. The real number Ta.f3[p\ formally defined by can be rigorously defined for any nonnegative 
function p such that ^ E Hp^^.{Q). 

2. IfV is constant, po — N/\Q\ is a stable local minimizer of \2.1]) . 

3. For V e L^/'^{Q) such that V - -J- / V is small enough (for the L^/'^{Q) norm), problem \2. if) 

1^1 JQ 

has a unique local minimizer in the neighborhood of po . 



4. Assume that V e LP{Q) with p > 3/2. When 

3/2 

An 




with 



the ground state energy In equals —00. 

It is easy to obtain a numerical value of Aq {Aq ~ 9.5785), hence of A^a,/3 for all {a, (3). The results 
are displayed on Fig One can see that the critical values Na^js are not very large. In particular, 

5 ± \/5 

Na^p ~ 4.636 for the values recommended in jg], namely = — - — . 
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Figure 1: Plot of the function a. i— > A^a,5/3-a for a. G [0,5/3]. 



Similar results can be obtained for models in which the electronic interaction is taken into account, such as 
the ones used in atomic-scale simulation of materials. Let us notably consider the minimization problem 



7;v = inf 



where p„ € L^^^{M.^) is a given nonnegative Q-periodic density such that / pn = N {pn represents the 

JQ 

density of smeared nuclear charges) , and where J and Ey^c are the Coulomb energy functional and some 
exchange-correlation energy functional, respectively. Recall that 

J[P - Pn] = / (P - Pn)W 
jQ 



with W denoting the unique solution in H^^^{Q) of 



-AW = 47r(p-p„), 



(2.3) 



For simphcity, we consider the case of the so-called Xa exchange-correlation functional 

JQ 



3 / 3\l/3 



where Cxc = is a positive constant, but similar results can be obtained for more complicated 

functionals. 
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Theorem 2.2 Let us consider two positive real numbers a and (3 such that a + /3 = 5/3. Then, 



1. If Pn — N/\Q\ (Jellium background) with N/\Q\ > pi^t — — I , then po = N/\Q\ is a stable local 

. L-TF , 



mimmizer o 



For Pn Q-periodic, such that / pn = N and close enough, for the L^^^{Q) norm, to some constant 

■Iq 

density po with po > pinf, problem \2.^) has a unique local minimizer in the neighborhood of pQ. 



3. When N > N^^p 



2 Ctf (gl^ - l) 



Cx 



3/2 

, the ground state energy In equals -co. 



Note that the constant pmf = I — — I — 0.102 is not very high (p = 0.110 for Al and p = 0.309 for Fe). 

\Ctf / 

3 Concluding remarks 

Density-independent orbital-free kinetic energy functional of the form lll.2|l are all nonlinearly unstable: 

5 ± \/5 

when the number N of electrons per unit cell exceeds a few units {N > 5 for a, P — — - — ), the kinetic 
energy goes to minus infinity when the density concentrates in some point of the unit cell. As proved in 
Theorem 1221 the Coulomb repulsion (which tends to prevent the density from concentrating) is not able 
to stabilize the model. For large, inhomogeneous systems simulated on fine grids, one can therefore fear 
that the numerical solution obtained with such models will be meaningless. 

This serious drawback is an additional motivation for constructing more elaborate functional such as the 
so-called density- dependent orbital-free functional introduced in 0. The mathematical analysis, as well 
as the numerical simulation of the latter models, are more difficult. Hopefully, this will be the matter of 
a future work. 



4 Mathematical proofs 

Let us begin this section by a formal calculation. In the sequel, the periodic lattice associated to the cell 
Q is denoted by TZ, and its dual lattice (see e.g. ^) by TZ* . If kp is a positive real number, and if / and 
g are two Q-periodic functions, one has 



Q 



f{x) I Wa.,i3{kF,x-y)g{y)dy] dx = f{x) > / Wa.^pikp, x - y - k)g{y + k) dy ] dx 



l{y)dy) dx = f{x) / Wa.pikp,- 

/ /(^) X! / ^aAkF,x~y-k)g{y)dy\ dx 

/ / f{x)g{y) \ ^Wa,f3{kF,x-y-k) \ dxdy. 
''Q-'Q \ken J 
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As TZ is symmetric with respect to the origin and as Wa,f3{kp, x) — Wa,p{kF, —x), it follows in particular 
that 



Wa,ii{kF,x-y)g{y)dy] dx ^ I g{x) I Wa.^pikp.x ~ y) f{y) dy ] dx 



and therefore that Ta^plp] — Tp^a[p\- In addition, using Poisson formula, one obtains 



^Wa,p{kF,x-y-k)^ ^ ■Wa,i3{kF,q)( 

ken qew 



i{x-y)-q 



9a/3 \Q\ 



^ ^ g(^ 

\Q\ ^^\2kF 



i{x-y)-q 



then 

fix) 



Wa.pikp.x - y) g[y)dy] dx 



9a/3 IQI K^kFj 



Cq{f)Cq{g) 



where {cq{h))q^Ti* denote the Fourier coefficients of the Q-periodic function h, namely 

Cq{h) = [ h{x) e-"'-'' dx. 



Proof of Theorem l2.H In view of the preceding remark, we can consider that ljl.2|l is a formal notation 
for 



JQ ^ JQ 



G 



TF 



9a/3 |g| 



q^n* 



2fco[p] 



(4.1) 



Note that the latter expression of the nonlocal term is the one which is actually used in numerical 

simulations (see e.g. [BJ). Now, it is easy to see that Tq_;3[p] is well defined by Ij4.1|l for any nonnegative 

5 

function p such that e H^^j.{Q) as soon as a and (3 are positive real numbers such that a -\- (3 ~ -. 

Indeed, when ^ € H^^^{Q), both p^^^ and ly^P are in L^{Q). Besides, the sum over the dual lattice is 
normally convergent. This can be proved by remarking that, on the one hand, G is a bounded function, 
and that, on the other hand, and p^ are in L'^{Q) for -y/p e H^^^{Q) and for a and /3 are in [0,5/3] 
(these are consequences of Sobolev inequalities QDl)- Thus, using Cauchy-Schwarz inequality and Parseval 
relation. 



- y 



G 



2fco[p] 



Cg (P") Cq [p^) < (^sup \G\^ 1^ iP"^ c« (P^^ 



qen* 



^ (»pi«i)^(EiM.-)r) 



1/2 / ^ 1/2 

2 



yqeTl* 



sup|G| ||p"|U.(Q)||p'^|U.(Q) <+oo. 
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This concludes the proof of the first statement of Theorem 12. II Let us now prove the second statement. 
For this purpose we introduce the problem 

mI^Ek[^]+ I^V^\ (/.€i/i„(Q), ^<^2^iv|, (4.2) 



where 



It is easy to check that the infimum of 114.211 is equal to In and that p is a local minimum of 112. l|l if and 
only if = is a local minimum of H4.2II . The Euler-Lagrange equation associated with problem (I4.2|l 
reads 

_iA0+^CTF|</'|'/'0 + y./) + i[0] =/i0 (4.3) 
where L[(j)] denotes the continuous linear form on Hp^j.{Q) defined by 

L[0]./, = Ctf^ 1^ ^^(2^^) ("c, {W^-'<kh) c, + /3c, (I0P") c, ilcbf^-'cbh)) , 

(4.4) 

and where ^ is a Lagrange multiplier associated with the constraint / (p^ = N. Let us denote by 



00 = \/Po- hand, (0o") = {^0^^ ~ ^ 9 7^ 0, and as on the other hand, 

G {0) = 0, one has L[(j)o] = 0. Therefore, (po solves the Euler-Lagrange equation 114. 3|l for V equal to the 
constant Vq, with fi = Ipq^'^ + Vq; if V is constant, (^0 is thus a critical point of 114. 2|l . In order to complete 
the proof of the second statement of Theorem 12 .11 it is sufficient to show that the continuous symmetric 
bilinear form 

B[cl)o,fio]{hi,h2) = l- f V/ii • V/i2 + ?Ctf / 0o^'/ii/i2 + 77^CTF/C[(/.o](/ii,/i2)-Mo / hih2, (4.5) 

is positive definite on the tangent subspace |/i G H^^^{Q), J (poh = oj-. In the above expression ^0 = 
Ipq^"^ and IC[(j>o] denotes the bilinear form defined by 

4a/3 ,2a+2/3-2 ^ / k 



4a/3 ^4/3 
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The latter equality has been obtained using that for all q ^ 0, Cq (^q") — (^'Po^^ ~ that G(0) = 0, 
and that a + /? = 5/3. A simple calculation then leads to 



B[(f>o,fio]{h, h) 



20 

9|Q| 

2 



qen* 



G 



\q\ 



3|Q| 

2 



3|QI 



G 



91 



G 



2fco[po] 
9 



|g| 
2fco[po 

1 - 3 



C, (/l)Cq (/l) 



2fco[po] 



Cg (/i)c, (h) 



3151 Coi^D'E/l.s^j 

where the function F is defined by 



1 1-772 



477 



log 



1 + r/ 



1-77 



Not surprisingly, one recovers the function ^(77) arising in Lindhard theory flj. As Firj) > 1 for all 77 > 0, 
one obtains 

B[^o,^^o]ih,h) > ^ (fcobo])' / (4.6) 

JQ 

which completes the proof of the second statement. 
Let us now consider the function 

: (H^^AQ) X R) X L3/2(q) ^ (i/p-i.(Q) X M) 

The function T is of class G^ and satisfies, for any constant Vo, -^((00,^0 + Vo) , Vb) — 0. Besides, the 
partial derivative of F with respect to the pair (</>, at the point {{(j)o, Mo + Vb) , Vb) is given by 



((0o.Po+yo),Vo) 



= (-00,0), 

((0o,Mo + Vo).Vo) 

where B[0o,Mo] denotes the bilinear form defined by 114. Next, it is possible to improve Il4.fi|l . by 
showing that there exists some constant 7 > such that ^(77) > 1 + 777^ for all 77 > (actually, one can 
use 7 = 1/4). Hence, we have 

B[cj,o,fio]{h,h)>l{ko[po]f f h^ + ll |V/7p. 
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This shows that i3[(/)o,/io] is coercive, and we may thus apply Lax-Milgram theorem 0j, proving that 
-^-r- — r- is invertible. The third statement of Theorem 12 . II then follows from the implicit 

function theorem |3]. 

The fourth statement can be estabHshed by a scaHng argument. We choose the coordinate axes in such 
a way that B{Q,e) ^ {x e R^, \x\ < e} C Q for some e > 0, and we consider a density pi € 

supported in B{0, e) such that pi > 0, / pi ~ N. We then consider the family of trial densities (po 
defined by 

Vcr > 1, Va; G Q, paix) = pi{ax). 
It is clear that for all cr > 1, belongs to the minimization set 



|p>0, ^ei/i„(Q), J^P = ^] 



One has 



qeK' 



2ko[p] 



c,{p^)cg (p^) 



As 



(p") = o- 



, and c„ 



we obtain 



Ta,l3[Po 



|VVPiI'+Ctf / P^ + C^TF 



R3 



9a/5 cr3|Q| 



2fco[p] 



Using that G is bounded and that lim G{ri) = — -, we have 



1 



^'1^1 ,S.^V2fco[p] 



1 



p? ('z)Pi (g) 



5 (27r) 



P? (?) Pi (5') 



5/3 
Pi • 



As TV > Na.,i3, and as C^(R3) is dense in H^{E.^), it is possible to find a function (f) € C^(M3) such that 



= 1 and 



|V0|^ 



10/3 



<2C.Ff^-l]^^/3 



(4.7) 
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For some t large enough, the function 

is supported in the set B{0,e) introduced above and the function pi{x) — N(f>T-(xY is such that pi G 

C^{R^), Supppi C B(0,e), pi > and / pi=N. In addition, 

Je3 



7 = 



One therefore has, 



IV^PTI' + Ctf 

Ta,t3[Pa 



9a/3 



Pil'/' >0. 



(7— > + CX5 



-7(7 



Besides, from Holder inequality @], 



/ 



<IIp.IIl.' \\v\\l. 



where p' = ( 1 - - ) < 3. As UpctILp' = o"^ IIpiUlp' = o{<7'^), we finally conclude that 



Ta.p[pa] + I PaV 
JQ 



a — *+oc 



-7Cr 



and therefore that In = — oo. 



Remark 4.1 Let us point out that one can carry out the same analysis with the density 

= |g| + Nccr'^piiax), (4.8) 

with Na,i3 < Nc < N instead of the above pa- This is physically more satisfactory since the densities 
defined by i4.8\) are uniformly bounded away from zero. 

Proof of Theorem 12.21 We use the same strategy as in the proof of Theorem 12. 11 and thus define the 
following minimization problem: 

inf|i?K[0] + ^J[</''-p„]+i?xc['^'], 0ei/p^,(Q), ^02^7v|. (4.9) 

A function (j) is a solution of H4.9II if and only if p = is a solution of H2.2II . Let us write down the 
corresponding Euler-Lagrange equation: 

_iA0 + ^Ctf|0|'/'</' + m +Wct>~ ^ae|0p/'0 = /i<^, (4.10) 

where the electrostatic potential W is defined by Ij2.3|l . and the linear form L[(j)] by Ij4.4ll . The constant p 
is the Lagrange multipHer associated to the charge constraint. As pointed out in the proof of Theorem l2.lL 
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if we define (fiQ — ^/po, we have L[0o] — 0. In addition, W is then a solution of —AW — 0, and is thus 
identically zero in view of its periodicity and of the normalization condition in 112.^^11 . This shows that 0o 
is a solution of Ij4.1flll with 

Hence, (j>o is a criticial point of the energy. We need now to show that it is a local minimizer. In order 
to do so, we show that the bilinear form 

z Jq ^ Jq yo^P 

20 f f 

+2(j)oD{hi,h2) ~ C^c^ <pl'^hih2- ^lo hih2, (4.12) 
y Jq Jq 

is positive definite on the tangent subspace |/i £ i/p(,j.(Q), J (j>oh — o|. In the above expression fiQ is 
defined by Ij4.11ll . /C[(/)o] denotes the bilinear form defined by 



and D is defined by 

D{hi,h2)^ [ Wih2, with - AT^i = 47r/ii, VFi e 7?ji„.(g), / VFi = 0. 

Jq 

We now point out that actually, Wi may be defined by its Fourier coefficients through |(7pCg(PVi) 
Aircq^hi) and co(Wi) = 0, so that 



47r Cq{hl)Cq{h2) 



Hence, carrying out the same computation as in the proof of Theorem 12. II we have 



qen-\{a} 1^' 

C 

Now, we know that F{ri) > 1, which, with the help of pj^'^ > , implies that 

Ctf 

for some positive constant S independent of h. This proves that po is a local minimizer of H2.2|l . 
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We now prove the second statement of Theorem [221 For this purpose, we introduce the function 



where W is here again defined by i2.l^ . where p = The function G is of class (all terms are clearly 
except W4>, but this one may be shown to have the desired regularity with the help of standard elliptic 
estimates). In addition, 5((0o, Mo), Po) = 0, and the partial derivatives of G at this point read 



dG 

del) 



((0O,A'o),Po) 

dG 



(h) = {B[cj)o,fio]{h,-), / M 



((0o,/io),Po) 



where i3[0O:Mo] is defined by II4.12|I . Here again, using the fact that F{r/) > 1 + 777^ for some positive 
1/3 

constant 7 and that pg > one easily shows that 



B[^o,^lo]ih,h)>5 



((0o,Mo),Po) 



n 9G 

for some constant S > 0. Hence, one may apply Lax-Milgram theorem H to prove that — 

is invertible. The second statement of Theorem 12.11 then follows from the implicit function theorem 
We now turn to the third statement of Theorem [221 We use test functions of the form 

We can carry out the same computation for the kinetic energy, showing here again that, if (pi satisfies 
H4.7|l . choosing pi{x) = N\(j)i{x)\'^ leads to 



(T — > + 00 



with 



7 = 



G 



TF 



9af3 



1 



|pi|'/^>0. 



We therefore only need to check that the remaining terms of the energy have a scaling of lower order as 
a goes to infinity. First, we have 



We then compute the Coulomb term, using its Fourier expression: 

47r 



\Q\ ^ \q? ~ \Q\ ^ 



\Cq{Pa)? + \Cq{pn)? 



13 



Pointing out, as in the proof of Theorem 12 .11 that Cq{pa) = Pi ) , 



we thus have 





) 



2 



J{Pa - Pn) < 



E 



+ c, 



qeK-\{o} 



where C is a constant depending only on p„. The sum is, up to a factor a, a Riemann sum, and we thus 



This allows to conclude that both the exchange term and the electrostatic term have a scaling of order 
strictly lower than cr^ as a goes to infinity. We thus come to the same conclusion as in Theorem 12. II 
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